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O ■ The matrix models which are conjectured to compute the circle Wilson loop and its 

correlator with chiral primary operators are mapped onto normal matrix models. A fermion 

p • droplet picture analogous to the well-known one for chiral primary operators is shown to 

pH I emerge in the large N limit. Several examples are computed. We find an interesting selection 

rule for the correlator of a single trace Wilson loop with a chiral primary operator. It can 

k> ! be non-zero only if the chiral primary is in a representation with a single hook. We show 

H ! that the expectation value of the Wilson loop in a large representation labelled by a Young 

diagram with a single row has a first order phase transition between a regime where it is 

identical to a large column representation and a regime where it is a large wrapping number 

single trace Wilson loop. 
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1 Introduction 

The basic gauge invariant operator in Yang-Mills theory is the Wilson loop. It is used as a 
diagnostic of the dynamical behavior of the gauge theory. The area law of the Wilson loop 
as a signal of confinement is a well known example. It has also provided an interesting tool 
in the AdS/CFT duality where it is the Yang-Mills operator which is the most direct probe 
of fundamental strings [1, 2]. 

Recently, highly symmetric Wilson loops have attracted considerable attention due the 
intriguing possibility that their AdS/CFT duals are described by symmetric 3-branes and 5- 
branes [3]-[10]. This is in analogy with large representation chiral primary operators of A/" = 4 
supersymmetric Yang-Mills theory which are interpreted as giant gravitons [11, 12, 13]. In 
the latter case, there is a beautiful picture of "bubbling geometries" which relates geometrical 
data of the space-time on the string theory side of AdS/CFT to the representations of SU(N) 
which are used to characterize chiral primary operators [14]. 

In that picture, the set of |-BPS chiral primary operators, the bosonic part of whose sym- 
metry group is 5*0(4) x 5*0(4) x R^, are Tr^(z(a;)) where R is an irreducible representation 
of the GL(N,C) Lie algebra. An operator has U(l) R-charge J if the Young diagram corre- 
sponding to R has J boxes. Then, there is a 1-1 correspondence between Young diagrams 
and Slater determinants, which are the wave-functions for quantum states of A^ free fermions 
in a harmonic potential. Due to the Pauli exclusion principle, the fermions occupy distinct 
states and in the appropriate large N limit, their semi-classical trajectories coalesce to form 
an incompressible droplet in the single- fermion phase space [15]. An analogy is then drawn 
between the shape of this droplet and the initial data that is necessary to specify |-BPS, 
5*0(4) X 5*0(4) X /?^-symmetric solutions of IIB supergravity, which are interpreted as the 
background space-times in which the giant graviton is embedded [14]. 

In this paper, we shall further develop the incompressible droplet model so that it can 
be used to compute the circle Wilson loop, multi-wrapped and multi-trace circle loops and 
correlation functions between circle loops and chiral primary operators. We do this by map- 
ping the Hermitian matrix model that is conjectured to describe the circular loop and the 
complex matrix model which describes the correlation function of a circular loop and a chiral 
primary operator onto normal matrix models. A normal matrix model is a model for complex 
eigenvalues. In the large N limit, the eigenvalue density forms a droplet in the complex plane. 
We argue that this droplet is identical to the one which characterizes chiral primary opera- 
tors and which describes the states of matrix quantum mechanics [15]. Our normal matrix 
model integrals simply compute overlaps between wave- functions of the matrix oscillator in 
a coherent state basis. 

One of the results that we shall find is that the operator product expansion of the funda- 
mental representation Wilson loop for a small circle W^q [circle] contains the chiral primary 



operator Ti jiz{x) only when i? is a representation of the gauge group with a single hook, 



for example R = l ■ Recall that, in the AdS/CFT correspondence, the Wilson loop is the 
source of a string world-sheet and the chiral primary operator is a graviton, or giant graviton 
if R is big. Our result implies that the integral of the vertex operator that couples a gravi- 
ton or giant graviton to the string world-sheet vanishes unless the graviton has a particular 
property. At present we do not understand the origin of this selection rule, but we do observe 
that, as a result of it, the Wilson loop operator acts as a probe of this aspect of the structure 
of a graviton. A single hook giant graviton is a single spherical 3-brane extended in either 
the AdS^ or S^ direction [13]. 

More generally, we find that the multi-trace circle Wilson loop operator (W^p [circle]) 
contains giant gravitons with up to K hooks. In the matrix model, this is an exact statement, 
holding for all values of A^, other parameters and representations R, though the reasoning 
that relates it to an exact identity for the Wilson loop is most reliable in the large A^ limit. 
It is a generalization of the same known result for the matrix elements of single trace chiral 
primary operators and giant gravitons, < Ttz'^{x)Tyii{z{0)) > is non-zero only if R has a 
single hook [16]. 

In the literature, there are two different descriptions of the Wilson loop which sources a 
D3-brane, in ref. [3] the claim was that the large wrapping number single trace Wilson loop 
whereas in [6] the loop was one in a representation with a long single row Young diagram. 
We will show that, if A is large enough, and in the large A^ limit, these two descriptions are 
in fact identical. The dominant contribution in an expansion of the character for the long 
row representation in symmetric polynomials turns out to be the large winding single trace. 
We shall also observe that, if A is lowered beyond a certain critical value ~ 5.5, there is a 
first order phase transition to a regime where the Wilson loop has free energy resembling the 
D5-brane as described in ref. [5, 6]. 

The loop operator of most interest in the AdS/CFT correspondence [1, 2], 



W[C] = ^TrPexp 



dr [tA^{x{T))x^{T) + 0(x(r)) ■ ^|x(r)| 



fl.ll 



is a measure of the holonomy of a heavy W-boson which is created by a symmetry break- 
ing condensate of the scalar fields in A/" = 4 supersymmetric Yang-Mills theory with 6'^ = 
(0*)/|(0)|. In IIB super-string theory, this operator provides a source for fundamental open 
strings with world-sheets bounded by the contour C, itself located on the boundary of 
AdS^ X S^. In (1.1) we have written the loop appropriate to Euclidean space which will 
be of most interest to us in the following. 

When C is a straight line or a circle, W[C] is a |-BPS operator. It is therefore con- 
jectured to be protected by supersymmetry. When C is a straight line, it is believed that 
(W^ [straight line]) = 1. This is supported by perturbative computations up to a few orders 
as well as the strong coupling limit computed using AdS/CFT. This is also the case for a 



multiply wound loop or a product of Wilson loop operators for any array of parallel straight 
lines. 

For the circular loop in Euclidean space, the expectation value at large N is thought to 
be given by the large N limit of the Gaussian matrix integral [17] 



2Nr 



JdM iTi e^) e-^TrM^ 2 

(W^[circle]) = ^^ > ^lA^fX) as iV ^ 00 , A fixed (1.2) 

where A = Qym^ i^ ^^e 't Hooft coupling. When computed in an expansion at weak coupling, 
at small A, this ratio is equal to the sum of ladder Feynman diagrams. It has been shown 
explicitly that the first few orders in diagrams which have not been included in the sum, 
those with internal loops, cancel identically [17, 18, 19]. It is conjectured, but not proven, 
that such corrections cancel to all orders and that the sum of ladder diagrams is the entire 
perturbative contribution. It is possible that instantons give a non-perturbative contribution 
[20, 21], so the sum of Feynman diagrams might not give the exact amplitude. However, in 
the infinite iV, 't Hooft limit, it is plausible that instantons are suppressed and the sum of 
the Feynman diagrams that is summarized in (1.2) is indeed exact. 

As further evidence for the consistency of this picture, when the limit of large A is taken 
on the right-hand-side of eqn. (1.2) it agrees with the AdS/CFT computation originally given 
in ref. [2], which should be valid in that limit. The circle can be obtained from a straight 
line by a conformal transformation. The fact that the circle and line do not have identical 
expectation values has been attributed to an anomaly [22]. In ref. [22], they also give a 
very nice argument relating, at large A, the -^ expansion of (1.2) to the moduli of Riemann 
surfaces in the topological expansion in string theory. 

As well as the fundamental representation loop in (1.2), the sum of ladder diagrams for 
higher representation loops is given by simply inserting the higher representation loop into the 
matrix integral. A similar equation applies to multiply wound Wilson loops or correlations 
of coincident loops 

(ri'*',i-iei>^ T^fMj^ <"^' 

Here, we emphasize, that this matrix model summarizes what is believed to be the sum 
of all Feynman diagrams. The complete set of |-BPS loops, the bosonic part of whose 
symmetry group is SL{2, R) x 5*0(3) x 5*0(5) is specified by considering loops in all irreducible 
representations of the SU{N) gauge group. The Schur polynomial formula can be used to 
express the expectation value of an arbitrary representation in terms of a sum over products 
of traces of multiply wound loops, the terms in which appear in (1.3). 

When probed from a distance much greater than its radius, r, a circle Wilson loop looks 
like a combination of local operators 

iy[circle] = (l^[circle]) ■ ^^a ■ (27rr)^OA(0) 

A 



where the center of the loop is near x = 0. We are ignoring the shght anisotropy which 
should occur in this formula. Since A/" = 4 super Yang-Mills theory is a conformal field 
theory, operators can be classified according to their conformal dimensions, A, along with 
other quantum numbers. The series contains both primary operators and descendants. The 
coefficients .^a of primary operators can be deduced from the asymptotics of the two-point 
function of the loop operator and the primary field, 

{W[ciTde]OAix)) ^ ^Oa{0)Oa{x))U ■ (27rr)^ + higher powers of l/x (1.4) 

( W [circle] ) 

Since the loop contains descendants as well as primary operators, the correlator with the 
primary operator is given by the leading asymptotic behavior at large |a;|. The descendants 
would appear in higher powers of l/\x\. 

For example, consider the chiral primary operator 

C*^*2---*'=(a;) = Tr ((l)'^{x)(f)''^{x) . . . 0*'=(a;)) symmetric traceless 

This is a |-BPS operator which has an exact conformal dimension A = k protected by 
supersymmetry. Two- and three-point functions of chiral primary operators are independent 
of the Yang-Mills coupling constant. We are most interested in chiral primary operators of a 
particular kind: if we take 6 = (1,0,...) in (1.1) and consider 

Oj{x) = Tr{z-^{x)) , z{x) =(l)\x)+i(l)'^{x) (1.5) 

which has exact conformal dimension A = J. 

The correlation function of this operator with a circle Wilson loop located at the origin 

has the form 

{W[ciTc\e]Oj{x)) ( 2nr -^-^ 



-] • O as |a;| ^ oo (1.6) 



(W^[circle]) ■ |0j|2 \4:Tx'^x'^ , 
The normalization of the chiral primary operator is defined by its 2-point function, 

{0.j{x)Oj.m = \Oj\Hjr (^)' (1-7) 

and can be computed using the complex matrix integral [16] 

/c^^^e-^^---Trz^Trz^' 
Jd^ze A 

The constants ^j in (1.6) were computed in the planar, small A limit using perturbation 
theory and in large A limit using the AdS/CFT correspondence in ref. [23]. The sum of 
planar ladder diagrams contributing to (1.6) was found in ref. [24]: 



IvOA/M) 
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This was conjectured to be an exact result. The sum of planar ladders which is summarized 
in (1.9), when extrapolated to large A, agrees with the AdS/CFT computation of the same 
coefficients ^j in ref. [23]. Leading corrections from other diagrams which are not ladders 
were also computed and shown to vanish in ref. [24], supporting the conjecture that they 
give the exact result. This has further been supported by comparison with large quantum 
number limits of strings in ref. [25]. 

The sum of planar ladder diagrams which contributes to (1.9) is also given by the large A^ 
limit of a matrix integral. As well as the single trace (1.5) we could consider a more general 
set of chiral primaries Tr^z(a;) where the trace is taken in an irreducible representation.^ 
In general, we could consider a Wilson loop in a representation R. Its expansion in local 
operators could contain a chiral primary in a representation R'. The coefficient C,rr' in this 
expansion can also be expressed as a matrix integral, 

^ /^^^dI^Tr^(e^(-+-"0 Tr^,(.) e-^^--- 
^^''' Id^^ 3^Tr« (e^(--)) e-^-- ■ |0^f ^''''^ 



dimR 

On the string theory side, chiral primary operators correspond to supergravitons and their 
coupling to loop operators can be extracted from asymptotics of the connected correlator of 
two far separated loop operators, which are sources for fundamental string world-sheets sitting 
at the boundary of AdS^ x S^. For small representations, they should give information about 
the coupling of gravitons to the fundamental string world-sheets that are sourced by the loop 
operators. For bigger representations they would yield information about the gravitational 
interactions of the more exotic objects, D3-branes and D5-branes, that are sourced by the 
loops. When R and R' are both large, (1.10) gives information about the coupling of giant 
strings to giant gravitons. 

In this Paper, we will take a closer look at some of the properties of the matrix integrals 
(1.2) and (1.10). In order to study them at large A^, we will rewrite them as models for 
complex eigenvalues, which are called Normal matrix models. 



^For an irreducible representation R corresponding to a given Young Diagram with J boxes, we begin 
with 



Z2 «3 «4 ^.^ %\ 



and consider hncar combinations of permutations of the upper indices until they transform in the represen- 
tation of the symmetric group that is associated with the Young Diagram. We then contract the indices 
to get a linear combination of single and multi-trace operators. For example, when J = 2, there are two 
representations 

z]zl -^ 2 V^J'^^ + ^i^O ^ 2'^^^^ ^ -TizTrz = TrQ-|Z 
z'jzf -> - i^z'jzf - z^^zlj -^ -Tiz^ - -TrzTrz = Tr^z 



1.1 Eigenvalue models 

In some cases, following the technique of Ginibre [26] , an integral over complex matrices can 
be written as an integral over complex eigenvalues where the matrices are diagonal and, with 
the suitable measure, we are left with integrating over the diagonal elements only. The most 
straightforward case is when the objects inside traces in the integrand are holomorphic in the 
complex matrix. An example is the integral which has to be done to find the normalization 
of the chiral primary operators 

Jd^ze-'^^^--TT^-zTT^,z 

The matrices can be made upper-triangular by a unitary transformation. The upper-triangular 
property (% = Hi > j) is preserved by products of z's and the trace of a product of the z's 
depends only on diagonal matrix elements, Tr{z'^) = J2i ^u- Iii the following, we will use one 
index to denote diagonal components, Zi = zu. The conjugate matrix, z, is lower-triangular 
and similar arguments apply. 

The matrix integral can thus be written as an integral over diagonal components times 
the appropriate Jacobian times an integral over the components in the upper triangle and an 
integration over unitary matrices. The integration of the latter two decouples and is common 
between the numerator and denominator in (1.11). They cancel, leaving the integral over 
diagonal components with the Jacobian. The resulting measure is 

N 

fl.l2) 



N-j 



When the matrix is triangular, the trace of a matrix in an irreducible representation R, Tr^z 
is a function only of its diagonal components and is given by the Schur polynomial 

T, , det(j-^"-^) det(.f-""'^) 

^'^"- det(.f-) - A(.) ^'-'^^ 

Here hi > h2 ■ ■ ■ > h^ are the lengths of the rows of the Young diagram associated with the 
representation R. 

It is now easy to see that that (1.11) can be written as an integral for diagonal elements 
of matrices, and then performed explicitly. 





d^z = 


i=l 


. ^iz)A{-z) 




with the vandermonde determinant 








A(^). 


Kj 


- ^.) = 


det 4~^ = 


det Zj 



{OrOr,) 



^ .j^j^ — j h'.+N—j 

/niI,^^..e-^N- |A(.)P '^^%^ '^%^ 
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E.,.,5.(-l)'^^^'^+^^^^ /nf=ici2^,e-^i-i^ (^f-'^^^^)!^^^^^^) 

where Sn is the set of all permutations of the integers (1, ..., A^) and we have decomposed 
the determinants into sums over permutations using the formula 

N 

detM= J2 i-lf"^''l[Mi^'^ 
In the large A^ limit, if we assume that hi/N -^ for all hi, the result of (1.14) reduces to 



'A^ 



E.'^^ 



IOrI"^ ^ ( ^ 1 d.s N ^ oo 



In the remainder of this paper, we are going to study matrix integrals such as (1-10) where 
making the matrices upper-triangular does not decouple the eigenvalues. What we will find 
is that, in many cases, we can nevertheless map it onto an eigenvalue integral over complex 
eigenvalues, that is, a Normal Matrix model. 

2 Normal Matrix Model 

A normal matrix model [27, 28] is a complex matrix model where the real and imaginary 
parts of the matrices are constrained to commute with each other. They can be simulta- 
neously diagonalized by conjugation with a unitary matrix, z = U ■ (iia,g{zi, ...,Z]^)U'^ and 
z = Udia.g{zi, ...,Zn)W where UW = I. Then, the integration of any function of the 
matrices which is invariant under {z.z) -^ {UzU\UzW) is taken as an integral over the 
diagonal elements Zi with measure given in (1.12) and the matrix Zij replaced by the diag- 
onal one diag(zi, ...Zn)- Normal matrix models have recently found a number of interesting 
applications [29]- [32]. They were discussed in the context of non-critical string theory in 
refs. [33]- [35]. 

The normal matrix model and the complex matrix model are equivalent if, as is the case 
in (1.11) solved above, aside from the Gaussian term, the quantities in the integrand are 
holomorphic in the matrices. However, in the case of the Wilson loop (1.2) or the correlator 
(1.10), the integrand is not holomorphic and it is less obvious that it can be written as an 
eigenvalue model. 



Nevertheless, there are circumstances where the normal matrix model is equivalent to 
the complex matrix model. In the remainder of this Paper, we will examine some of these 
examples. The interesting feature of a normal matrix model is that, in the large A^ limit, the 
distribution of eigenvalues has support on a two-dimensional subset of the complex plane. 
This will lead to a droplet model for the eigenvalues. 

2.1 Fundamental representation Wilson loop 

As an example, consider the expectation value of the Wilson loop in the fundamental rep- 
resentation which is given by the Hermitian matrix integral (1.2) of the Wilson loop in this 
representation. 



J[dM] e-^TrM^ ^Tr e 



,M 



We can easily write this as a complex matrix model. The approach (which we shall not 
follow) would be to make this into a complex matrix model by introducing the imaginary 
part as another Gaussian matrix Mj, so that, z = M + iMj and writing 

(Wa) = ^r (2.2) 

It is not obvious how to now write (2.2) as an integral over eigenvalues. Backtracking, we 
recall that we could also have written (2.1) as an integral over the eigenvalues of Hermitian 
matrices, 

.... . _ /n.c/m,e-^-?A^(m) (e^Q 
J Yiidrriie a ™»A^(m) 

To proceed, we realize that the integral that we must do in (2.3) is formally the expec- 
tation value of the operator e"^^ in the quantum state given by the wave- function iplrrii) = 
A[m)e~~ ^i "*' . This is just the ground state of a system of A^ fermions moving in a Harmonic 
potential well. (To see this, we must recall some properties of determinants of polynomials. 
This point is discussed in the Appendix at the end of this Section. See eq. (2.19).). There 
are other ways to present this integral. One of them is uses a holomorphic polarization of the 
oscillator phase space. In that presentation, the ground state wave-function of an N-fermion 
system is ip{z) = A{z)e~~ l^i^^^^ where z is a complex variable and the expectation value is 



lWid^Zie-^-^~'^'^\A{z)\^ ipTr (eVl("+^") 
i^M^) = e-V8A^ V y_ (2.4) 

A detailed derivation of a general integral formula of this kind, eq. (2.17) is given in the 
Appendix to this Section. 



This is now a normal matrix model, written as an integral over complex eigenvalues. 
Note two differences between the normal model (2.4) and the complex model (2.2): the 
factor e^^/^^ in front of (2.4) and the factor of 4^, rather than | in front of {z + 



V2' 



Z . 



In summary, we have found three identical presentations of the same matrix model: 
{Wa) = -^ l[dM]e~'^^'^'^Tie^ (hermitian model) 
= ^ /M^^] e-^T'-^^^Tre5(^+^) (complex model) 
= -^e-^ I [rf^z] e-^^''"^^Tre75(^+^) (normal model) (2.5) 

Z]^ J[2,2]=0 N 



where we have abbreviated the denominators as the Gaussian partition functions Zh-, Zq-, Zn 
for Hermitian, complex and normal matrix models, respectively. We can easily check that 
the order A term agrees in all three computations in (2.5): 

I , ^, I \ N^ A 

= l + ^(TrMV.n..an + --- = l + ^^- + --- = l + - + ---, 

]^ 1 A 2iV^ A 

= l + i]7^TrZZ)_piex + --- = l + ^^^ + --- = l + - + ---, (2.6) 

= l + ^(TrZZ)™,-A + ... = i + ^A^^:!±i^_A + ... = i + ^ + .... 

2N^ /normal g^ ^ 2N 2N 2 8N 8 

We also emphasize that the three expressions in (2.5) are exactly equal for all values of A 

and A^. 

2.2 Fermion droplet model 

Let us now consider the multiply wound circular Wilson loop with k windings. The discussion 
of the previous Section leads to the Normal matrix integral 

Zjv J[z,z]=0 N 

This can be written as an integral over diagonal matrices with 

(W"n ^) = ^^ e- w / J] d^z^ exp — - f^ ZiZi + ^ In \zi - Zj\^ + -=(2:1 + z^] 
^N J i y A .^-^ ^^. V2 

In the large A^ limit, the integration over the eigenvalues ^2, ..., -Zat is dominated by a saddle 
point where the eigenvalues satisfy the saddle point equation 



1 ^ 1 h 



10 



If we introduce the normalized eigenvalue density 



1 ^ r 

p{z) = j^—^Y.^\z-z,) , jdPzp{z) = l 

the saddle point equation 

2 A^ - 1 /■ ,2 , p{z') 11 , , „, 

-. = ^^y rf . — -, + -— - . G support of p (2.7) 

There are two different regimes that we can analyze: 

1. Small winding loop: limAr_^oo Ji = ^ 

In this regime, the last term in the saddle point equation (2.7) is irrelevant and the 
equation is solved by the constant density droplet 

p{z) = \^>^ V 2 (2.8) 

There are two ways to compute the expectation value of the multiply wound Wilson 
loop. We can first simply use the formula for evaluating the expectation value of the 
trace of matrices once the density is known: 

(iy(')) = f d'zp{z)e^^"^'^'^ = \ f^rdr ^ dO gv^^^-^^ 
J nX Jo Jo 

= 1 f\dt r dOe^'^"''''^ = f\tdt Io{ky/\t) = -^h{k^) (2.9) 
TT JO Jo Jo sjXk 

This is valid when k « N and it agrees with the Hermitian Gaussian matrix model 
(1.2). 2 

Alternatively, we can use the large N limit to integrate out the variables Z2.,....,zn to 
get an effective theory for zi interacting with the droplet. The result is an integral for 
Zi with effective action 

S'eflf = —^^i^i - (^ - 1) / d^z\n \zi - z\^p{z) + -j={zi + zi) 

This action has an electrostatic interpretation. The fist term is the potential due to a 
constant charge density on a large disc and it tends to attract the oppositely charged 



^The expectation value of the single trace Wilson loop is known explicitly at finite A^ [3]. A useful formula 
is a contour integral representation [36] 



<»o>=^/2S^^Mt^i •'■'"' 
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particle with position Zi to the center of the disc. The second is the repulsive two- 
dimensional coulomb interaction with the charge distribution, which, using (2.8), is also 
a charged disc centered on the origin. We expect some cancellation of the attraction in 
the first term and the repulsion in the second term, particularly in the interior of the 
discs. The third term is the contribution of an external electric field which pulls the zi 
to the right in the complex plane. Using the constant density droplet. 



d'wp{w)\0g\w-z\' = <(aI^I+^^^2 FI^V2 (2.11) 

InlzP \z\ > J^ 

\ II I I V ^ 

we arrive at 

^ rfzi^,-(iV-l)(lnf-l) + ^(^i + ^i) + O(iV0) |^ip<| 
"'' I '-f-ZrZ,-{N-l)\n\z,\^ + f^{z, + -z,) + 0{N') |^ip > | ^''''^ 

Since we only found the eigenvalue density to leading order in large iV, only the order 
A^ terms in this effective action should be trusted. In the region |zip > I, the effective 
action has a minimum at its minimal value, |zip = I, and a steep slope, of order N 
as \zi\ increases. We conclude that, in the large A^ limit, the eigenvalue is confined to 
the disc \zi\ < ^. This is just the same disc where the eigenvalue distribution (2.8) has 
support in the infinite N limit. Inside the disc, the terms of order A^ in the energy are 
constant. To that order, the potential there is flat and, to understand the behavior of 
the particle, it is necessary to understand the terms of order one. However, to correctly 
evaluate the order one terms in the action (some of which have been written in (2.12)) 
we would need to properly evaluate the 1/N corrections to the density. If we did so, we 
would find that the terms of order one are also constant and the potential inside the 
disc is actually flat. The average of the position of the charge at Zi over the degenerate 
saddle point, which is the disc, and of course gives precisely the same result as (2.9). 

2. Large winding loop: limAr_,oo j^ = ^, with ^ a constant. 

In ref. [3], it is argued that the loop W^ with large winding k is dual to a D3-brane 
whose world-volume is AdS2 x S"^ embedded in AdS^ and with k units of electric flux 
on the D3-brane. It is argued in [6] that the AdS2 x S"^ D3-brane corresponds to a 
Wilson loop in the symmetric representation W^n^.-.n. 

In the following, we will analyze the large winding number limit of a single trace Wilson 
loop. In particular, the relation between D3-brane and large winding loop is checked 
in the limit 

k.N^oo with K = fixed (2.13) 
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In this case, the saddle-point equation (2.7) still holds and again the solution at leading 
order in large A^ is the disc-shaped droplet (2.8). In the effective action (2.12), there 
are now terms of order A^ in the action inside the disc as well as outside and it is 
these terms which must be taken into account to get the leading order in large A^. The 
integral over zi has the form 






|2i|2> 



2^2i 



rf22ie"^l''l+'^'°l''l +"^72(^1+-"!) (2.14) 
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In the large A^ limit, the classical position of the charge zi is outside of the disc. The 
second term in (2.14) dominates and can be evaluated by saddle-point integration. The 
saddle point is 

zl = ^(K + VlT^) = ^e^-'-^'^ (2.15) 

The integral is 

{W^'^) ~ e2N{nVJ^+sint-^.) ^ ^ = t^ (2.16) 

This expression agrees, as it should, with one derived in ref. [3] using the Hermitian 
matrix model. It was also shown in ref. [3] that this agrees with the Born-Infeld action 
of the D3-brane. It was further argued there that, for a large winding number Wilson 
loop, the fundamental string world-sheet blows up into a D3-brane. 

In our matrix model picture, the D3-brane corresponds to an isolated eigenvalue residing 
outside of the droplet. This has a strong analogy with the bubbling geometry where a 
giant graviton, which is also a D3-brane is also described in this way [14]. 

3. intermediate case: k ~ A^" with < a < 1 

In this case, the electric field due to the Wilson loop term in the action is not strong 
enough to pull the charge out of the disc, but in the limit A^ — » oo it is pulled to the 
edge where zi = zi = J^. 

The essential difference that we have found between the small and large winding number 
Wilson loop is that the eigenvalue which interacts with the loop resides inside the droplet for 
small winding number and outside of the droplet for large winding number. 

Appendix A: An integral formula 

In the following, we shall prove the integral formula 



/ nil dx,e--'A{xf nil e^'"' _ / Uli d^z,e-\'^\'\A{z)\^ Uli e^^^'+^'^'^'' 
IIlf=idxie~-'A{x)^ ~ JUf=id''zie-\-^\'\A{z)\^ 
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(2.17) 



which relates the expectation value of a general exponential in the Hermitian matrix model 
to one in the normal matrix model. 

The matrix eigenvalue integral on the left-hand-side of eq. (2.17) is written as 






r ^ 
Jl[dx,^{x,,---,XNfl[{(^'''"') (2-18) 



where \l/(xi, ■ ■ ■ , xn) is the Slater determinant which is the wave-function for the A^ fermions 
in a harmonic potential 

^(a;i,---,X7v) = det (V^i_i(a;^)) 

l<t,j<N 

^x) = ^i i/n(x)e-^"' (2.19) 

Hn{x) = e^ {—dx)"e~^ is the Hermite polynomial. In (2.18), we used the relation 

l^<ij(x„...,x^r = -^e- ^l^ -? A(x)2 (2.20) 

where Z^ is the partition function of hermitian matrix model 

N o N 



Zh = Jfldx,e-^'A{x)^ = 2-^(27r)f f[ n\ (2.21) 



i=l n=l 



We can further rewrite (2.18) in the operator formalism using the N free oscillators ai,al 
{[ai,aj] = 6ij). Using the relation Xi = ^^ , (2.18) becomes 



/njIi^a^»e-"'A(x)'niIie^'"' _ 1 ^ 

lulidx.e-^^Aixy " m ,=1 

where |\1/) is the anti-symmetric A^-particle state 



^1 J|e5'=»('*'+'*I)|^) (2.22) 



1^) = |0) A|l) A---A|A^-1) 

= E (-l)'''''"l^(l)-l)®---®l^(^)-l) eHi®---®HN (2.23) 

and \n) = -r7('2^)"|0) is the n^^ excited state. Here Tij denotes the Fock space oii^^ oscillator 
and the operators ai,al act on ?ij. The state |\E') represents a Fermi sea of N fermions in 
a harmonic potential. To evaluate (2.22), we need to know the overlap (n|e^^|m). In order 
to rewrite the hermitian matrix model to normal matrix model, the key step is to insert the 
completeness relation of the coherent state 1^;) = e^" |0) 

1 = - [ d^ze-\'\'\z){z\ (2.24) 
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Then the matrix element {n\e^^\m) is written as 

(n|e ^ \m) = e * (n|ev^ e\/2 |m) 



-ife^ 1 /■ ,2 -UP/ I -^a, w I 4=at| V 

e -i — a ze ' ' (n|ev^ |2;)(2;|ev^ |m) 

TT J 



^ ^rf2^e-l^l'+^(^+^)-^'^"z"^ (2.25) 



Let us consider a term appearing in (2.18) of the form 

. N N N 

/ l[dx^^ixu■■■,XNyl[e''^''^ = (^iHe^'^l*) (2-26) 

Using the relation (2.25), this is written as 

T N 1 N 

— (M/ine'^^^l^) = ]^ E (-l)'^+^nKO-l|e''^'k«-l) 

1 A' <. i- 



E (-i)^+^n/^v-'"''^^^ 

— n /rf'^.e-l^'l'+^(^^+^')-^'?A(^)A(^) (2.27) 



N 

n 

where Zj\f is the partition function of normal matrix model 

N N 



Zn = Jfl dhie-\'^\'\A{z)\^ = TT^ n n! (2.28) 



i=l n=l 



From (2.20) and (2.27), we find the identity between the Gaussian hermitian matrix model 
and the Gaussian normal matrix model 



^ l[dx,e'^"A{xfY[e'^^^ = -^ / J] rfV-l^'l>(^)r H e^(^'+^')-*'» (2.29) 

We should stress that the relation (2.29) is exact. 

3 Antisymmetric representations 

We showed in the previous section that, when the winding number k of the loop was small, 
it was described by an eigenvalue moving in a flat potential in the interior of the droplet. 
When k was larger, of order A^, the particle moved outside of the droplet and is interpreted 
as a D3-brane. In both cases, the expectation value of the Wilson loop was given by a certain 
electrostatic interaction energy of the particle with the rest of the droplet. 

In the dual string theory, the small winding loop is well described by the boundary of 
fundamental string in AdS^ x S^. We will now consider the situation when large number of 
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such small winding loops are put on top of each other. For definiteness, let us consider the 
loop (VFak) i'^ the K^^ anti-symmetric representation. In refs. [5, 6], it is argued that Waj^ 
is dual to a D5-brane of the shape AdS2 x S*^ with the K unit of electric flux on it. 

The K^^ antisymmetric representation is the Young diagram with one column of K boxes, 

Ak = -(^ boxes 
and we consider the expectation value of the Wilson loop in this representation, 

,,,, , J[-iA»-^""'^igTV.. (^") ,3., 



where dim^x = N\/K\{N — K)\. ^ This can be written as an integral over the eigenvalues 
of Hermitian matrices, 

{Wa^) = (-1)"^+' _g^ , ^ '- (3.2) 

JYlidrriie a ™»A2(m) 

Where we have used the Schur polynomial formula which gives us 

Tr^^ (diagM) = (-1)^+1 Y^ M,,M,,...M,^ 

ii<i2<---<iK 

(See Eqn. (1.14) for the example of -ft' = 2.) Note that, in Eq. (3.3), the combinatorial factor 
that we would find from choosing the first K eigenvalues cancels with the dimension of the 
representation. 

Once the Hermitian matrix integral is written in the form (3.2), we can apply the integral 
formula (2.17) to write the rewrite it as a normal matrix model:^ 

(Wa,) = (_i)^+ie-^A/8A. ^ J_ (3 3) 



^For a general representation R with Young tableaux (hi, • • • , h^), the dimension is 

hi — hj + j — i 



dimi? = Yl 



l<i<j<N -^ 

^Note that there is an interesting K ^ N — K duality in the integral (3.3). We find it by translating all 
of the eigenvalues mi ^> mi + \ r^h ^nd noting that the Vandermonde determinant is invariant under the 
simultaneous translation of all eigenvalues. Then, we transform mi -^ —mi to obtain 

In particular, (Wam) — ^^ ■ Since An is trivial as a representation of SU{N), {Wan) is equal to the circle 
Wilson loop in the non-interacting adjoint U{1) gauge theory. 
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To analyze this model, it is convenient to divide the eigenvalues into two sets, those that 
appear in the loop zi, ...,zk and those which do not, zk+i, ..., -Ztv- 

Let us first analyze the case where K remains finite as A^ ^ oo, i.e. limTv^oo ^ = 0- ^^ 
this case, the set oi N — K eigenvalues can be integrated out by solving for their saddle-point, 
which is again the circular droplet (2.8). The remaining K eigenvalues have an effective action 
which concentrates them within the circular droplet. Inside the droplet, the terms of order 
A^ in the effective action cancel, and to analyze the dynamics, one must take into account 
the 1/A^ corrections to the density, to get the effective action to order one. A guess would 
be that the potential for eigenvalues inside the drop is simply flat so they move freely there. 
Then, our result for infinite A^, finite K would be 

(Wa.) = i-lf^HW^r (3.4) 

We shall see that this is borne out by a careful analysis below. 

The effect of integrating out the N — K eigenvalues is summarized by the joint probability 
distribution of remaining K eigenvalues [26] 



Pn\.Zi, ■ ■ ■ , Zk) 



N\ I ntK+i d'^^ e-^ Et. i^^p n.<, k. - ^.f 



(^-^)' /nil.rf^^.e-^Et.N^n. 



i<j \^i 



When A^ -^ oo, we must study the convergence of the sum in (3.5). It is easy to see that the 
single particle distribution converges to a positive constant if \zi\ < J^ and zero if \zi\ is 
outside of this radius [26]. This means that the eigenvalues are confined to a droplet. Then, 
it can be seen that the multi-particle distribution converges to a constant for all eigenvalues 
inside the droplet. In this way, we see that the eigenvalues are confined to the droplet and 
have a flat potential when inside the droplet. 
In the large A^ limit with fixed K^ p becomes 




p{zu---,zk)={^) e'^«^i-l^»l det(e-^'^0^ ^ \^^\ < \h (3-6) 



It is easy to see that this produces (3.4). We will derive this result in a third way in the 
Appendix to this Section. 

In the case where K ^ oo and the ratio K/N remains finite as A^ — > oo, we again divide 
the eigenvalues into two subsets, so that one subset is those which couple directly to the 
Wilson loop and the other is those which do not. In this case, each subset contains a finite 
fraction of the eigenvalues. We introduce two densities, 

P(^) = 4 E ^(^ - ^^) ' P(^) = j;A^ E ^(^ - ^^) / d'zpiz) = 1 = 1 d'zpiz) (3.7) 

j=l i=K+l 
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and the saddle point equations are 

— 2; = -77 / S'z 7 H — — / c?^/ 7 H 1= — z G support of p 

A N } z-z' N J z-z' y/2N 

—2; = — — — / (fz + -r- / (fz z E support of p 

A N J z-z' NJ z-z' 

In the large A^ limit, these would apparently be solved by any two functions such that 

-p{z) + -^^p{^) = ^' ^^ 3.8 

They must add to form the circular droplet density in Eq. (2.8), but otherwise any split of 
the disc distribution into two functions of the appropriate normalization will do. This is a 
highly degenerate solution and, again, the degeneracy must be split by corrections of order 
1/A^. Here, it can be resolved by maximizing the integral 



E;^t = K J d^zp{z){z + z) 



In the electrostatic analogy, this term is the interaction of a subset of the charged particles 
with an external electric field. A fraction K/N of the particles interact with this field and the 
remaining fraction 1 — K/N do not. This polarizes the droplet by pushing the K eigenvalues 
whose density is p to the right of the droplet and the other N — K eigenvalues to the left 
(see figure lb). The result is 

P(^) = f ^ (v^- 1^1 I 0(^ + ^- V2Acos^) (3.9) 

P(^) = JfTJ^ix ® I V ^ " '"^'l ®^^^^ cose -z-z) (3.10) 

where 6 is the Heavyside function and the normalization of the densities in Eq. (3.7) deter- 
mines 9: 

1 K 

e- -sm2e = TT— 

2 N 
Then the expectation value of the Wilson loop is given by 

{Wa,) = e^v{^j d'zp{z){z+-z)\=e^^l^^^sin^d^ (3.11) 

which agrees with the result of Yamaguchi [5] who argued that it matches the free energy of 
the 5-brane computed using the Born-Infeld action and who derived the same result from a 
Hermitian matrix model. This agreement between hermitian model analysis and the normal 
model can easily be seen directly by projecting the circular droplet to the real axis of z. Then 



the eigenvalue density along the real axis is given by the semi-circle law, which of course is 
a eigenvalue distribution of the Gaussian hermitian matrix model. 

We summarize the results so far in Figure 1. For the small winding, small representation 
loop, the eigenvalue which couples to the loop lives in a fiat potential inside the disc. When 
the winding number gets large, it acts like an electric field which pulls the eigenvalue outside 
of the disc as depicted in fig. la. This is associated with the D3-brane. In the case where 
there is a large multiplicity of small winding number loops, the eigenvalues still reside in the 
disc, but the loop polarizes them as in fig. lb. One can think of the line which divides the 
two subsets of eigenvalues in this case as the analog of a hole in the eigenvalue sea. In the 
case of a chiral primary operator, such a hole is also interpreted as a giant graviton which is 
extended in the S^ direction. In the case of the Wilson loop, it seems to be more economical 
to interpret it as a 5-brane [6]. 





a) Large Winding Loop 



b) Antisymmetric Rep 



Figure 1: Eigenvalue distribution in the normal matrix representation of Wilson loops, a) A large 
winding loop W^''' corresponds to a single eigenvalue sitting at z = z* outside the droplet, b) 
For a Wilson loop in the K^^ antisymmetric representation Ax, the droplet is divided by a line 

Kez = y 2 cos 6 and K eigenvalues and N — K eigenvalues are distributed on the right and left of 
this line, respectively. 



3.1 Symmetric representation 

There seem to be two competing descriptions of the D3-brane. In ref. [3] the claim was 
that the large wrapping number single trace Wilson loop is dual to a D3-brane. In [6] the 
Wilson loop was one in a representation consisting of a long single row Young diagram. In 
the following we will show that there are in fact two regimes. If A is large enough, the two 
descriptions are in fact identical in the large N limit. If A is lowered beyond a certain critical 
value ~ 5.5, there is a first order phase transition to a regime where the Wilson loop behaves 
like a source of D5-brane as described in ref. [5, 6]. 

We can also consider the Wilson loop in a representation which has Young diagram a 
single row iinii...n with k boxes. The decomposition of the character in this representation 
into symmetric polynomials, for a diagonal matrix is 

K\{N -l)\ 



Tr 



M 



nin...n<: 



{N + K-l)\ 



E 



^mij+mij + .-.+mij. 



(3.12) 



ii<ii...<ik. 
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In order to use the formula (2.17), it is convenient to divide the sum in Eq. (3.12) into special 
cases. Here, we will compare the two extreme such cases, first the set of terms where all of 
the indices are the same. This corresponds to a single large winding Wilson loop. Second, 
we will consider the case where all indices are different. This gives a contribution identical 
to the antisymmetric representation Wilson loop. In both cases, we will consider the case 
where k is large. 

One would naively expect that the case with all indices different would dominate simply 
because it has a much higher multiplicity. In fact, if we were averaging the right-hand-side 
of (3.12) under an integral which was symmetric in the indices, the terms where all indices 
are different could be replaces by 



k\{N -k)\ 



miH-m2+...+mfc 



The factor comes from the multiplicity of these terms. 
The other term of interest would be replaced by 

and it occurs which a smaller factor. 

The expectation value of the first of these is identical to the expectation value of the 
Wilson loop in an anti-symmetric representation which we have already computed. Taking 
into account the multiplicity factor and using Eq. (3.11) we get 

where ^ — ^ sin2^ = -kjj. 

This should be compared with the exponential in N part of the expectation value of A^e^™^ 
which we can read off from Eq. (2.16), which gives 

kVx 



exp (A^ (2k^k'^ + 1 + 2 sinh"^ k)) , « = ^^ (3-14) 

To compare these, we choose ^ = \ and plot the difference between the exponents in 
Eqs. (3.13) and (3.14) as a function of A. The plot is depicted in figure 2. There we see 
a crossover between two regimes: when A < 5.5 the case (3.13) dominates whereas for 
A > 5.5 the case (3.14) dominates. We have not analyzed intermediate cases where sum but 
not all eignevalues are equal. We speculate that there is a phase transition between these 
two regimes. Seen from the point of view of the electrostatic analogy, in the first phase the 
eigenvalues which couple to the external electric field prefer to be independent. In the second 
phase they prefer to sit on top of each other as one macroscopically charged particle. This 
occurs when A is sufficiently large. 
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Figure 2: This is the graph of the differences in the exponents in (3.13) minus (3.14) (divided by 
N) in the region where they cross. There is a first order phase transition where A ~ 5.5. This figure 
is plotted for the case where k = N/2. 

Appendix B: Contour Integral Representation 

In this Appendix, we rederive (3.4) by directly integrating out the N — K eigenvalues using 
the joint probability distribution (3.5). 






K 72 



2_l1, /a 



N f^- 



TT 






AP 



(3.15) 



Here and the following section, we use the rescaled variables z'!^'^'" = \/^zf'^ to make the 



A "^i 



Gaussian part l^f'^'^P simple. 

As we will explain in the next subsection of this Appendix, we can integrate out the radial 
coordinate of Zi and the the angular integral reduces to a contour integral. 
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d'^Zi 
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-\^^\^+h^/w(^^+^') 
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K P-(^)jrPi 
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I<A -r-r 

^\J 2mz, 



1. /A 
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= TTr TT^esM TT i ^e2^^^> det — ^^ — ^-^ 3.16 

{N-K)\ l\f2m 1<^,<K[ ,^_,^^g ) 

We should stress that this is an exact expression for anti-symmetric representation at finite 
A^ and K. In the large N limit with finite K, the dominant contribution comes from the 
diagonal elements of the matrix inside the determinant. In this way, we see that the integral 
factorizes 






^-2^' - 



"* , 



- (^)''"n/ij-----'^0(VX))^(H^o)'>- (3,17) 

Although (3.16) is an exact formula, it seems difficult to take the large A^ and K limit with 
fixed K/N in this expression. 

Integrating Out the Radial Coordinate \z\ in (2.17) 

As promised above, let us consider the integration of the radial coordinates \zi\ of the normal 
matrix eigenvalues Zi. The remaining angular integral can be written as a contour integral. 
In (2.25), the matrix element of e v^ " was written as the integral over z in the coherent 
state \z). We write this as an integral over \z\ and arg(z) = 9 

\/n\m\ h Jo 2tc 

dt I ^e-*+^^(^+^"^)-^'='t^^"- (3.18) 



y/n\m\ Jo J 27riz 

where we introduced the variable t = |2;p and rewrote the ^-integral as a contour integral. 
After rescaling z —^t~2z, the t-integral can be performed easily 

Vn!m! -'o J 2'Kiz 

m\ f dz J^z-^k'^ z"^ 
eV2 4 





n! J 2m {z - ^)™+^ 

ev^'+* ^ 3.19 

2m z™+i ^ ^ 

In the last step, we have shifted z ^ z + -^. Using this expression, our key formula (2.27) 
is written as a contour integral 

1 ^ ^ 
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1 N 1 f ■y J- ^a \(j(a)-l 

CTjTgiijv a=l " ^a 

This almost looks like a unitary matrix model. For example, when 

ki = ■ ■ ■ = kx = ^, kx+i = ■ ■ ■ = kx = (3-21) 

(3.20) becomes a coupled system of U{K) and U{N — K) unitary matrix model 






j=i ^'- i=i-"^ ^" J 



^h 2ix -^ h 111 



X J] n (e'^" + ^ ~ e^'^0(e"*''' - e-^'^0 (3-22) 

i=l j = l V 2 

4 OPE of Wilson Loops and Chiral Primary Operators 

So far the normal matrix model representation of Wilson loop and the droplet picture which 
emerges at large A^ are interesting technical tools, but all of the computations that we have 
outlined could also have been done, in principle, using the Hermitian Gaussian matrix model. 
However, when we consider the interaction between Wilson loop and chiral primary operator, 
the normal matrix model representation provides us with new insight. 

As discussed in ref. [14], a 1/2-BPS chiral primary operator has a dual description in 
supergravity which is specified by the shape of a droplet. If we think of the chiral primary 
operator as a state, the characteristic of that state which is visible in geometry is the shape 
of the droplet that is associated with it. 

The disc-shaped droplet that we have discussed so far is the ground state. The state of IIB 
string theory corresponding to it is the ground state on the AdS^ x S^ background. Small 
chiral primaries correspond to the same background with some gravitons excited. Larger 
chiral primaries are giant gravitons or, if they are large enough, they modify the geometry 
and the state is string theory on a different background. In the large A^ limit these states are 
all characterized by droplets which, now are different from the ground state disc droplet. 

The vacuum expectation value of the Wilson loop operator measures the vacuum matrix 
element and gives the content of the unit operator in the loop. We can consider some other 
matrix elements as well, such as the matrix element between the vacuum and the state which 
is created by a chiral primary operator. This matrix element is of course given in the integral 
in eq. 1.10. In the language of eigenvalues, it measures the ability of the Wilson loop operator 
to modify the circular droplet which characterizes the ground state so that it fits with the 
droplet for an excited state, the chiral primary. 
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Let us consider the OPE between the Wilson loop and the chiral primary operator Ok = 
-jl=TiZ^. In the next Section, we shall prove that it is given by a normal matrix model 

The normalization of Ok is chosen so that the two-point function is normalized to unity in 
the large A^ limit. In the eigenvalue basis, (4.1) is written as 

(W^oO.) =2-|^/nrfV-l-nA(.)p-^ E e^v^(-+^^)4 (4.2) 

This is evaluated exactly using the joint probability distribution (3.5).^ 
As discussed in the previous section, we can integrate out the radial coordinate \zi\ and 
rewrite this OPE coefficient as a contour integral. However, to see the relation to the droplet 
picture, we find it useful to integrate out the angular part of Zi. After some computation, we 
find the exact finite A^ expression for the normal matrix integral (4.1) 



In the large A^ limit, the |z|-integral is localized to the saddle point \z\ = yN. Physically, 
this corresponds to the localization of the eigenvalue to the boundary of circular droplet. 
This is consistent with the picture that the chiral primary Ok with k <^ N corresponds to 
a graviton which in turn is represented by a small ripple of the Fermi surface. Finally, the 
OPE coefficient at large A^ is found to be 

i^nOk) = ^2~-^VkIk{VX) (4.4) 

which produces the expression in Eq. (1.9) which was originally derived in ref. [24]. 

In the derivation of (4.4) we have assumed that k <^ N. One interesting problem where 
this is not the case is the BMN limit k = J r^ yN . In this limit, the integral (4.3) behaves 
differently. To take the BMN limit, it is useful to rewrite the integral (4.3) as a sum of 
Laguerre polynomials 



'"'This is a more subtle application of the joint probability distribution since it probes its next-to-leading 
order behavior in 1/7V. The naive leading order in large N , which would be the factorized product of two 
traces, vanishes after the integral over phases of z. This means that we must take into account the next- 
to-leading order. We do this by integrating the exact formula to get (4.4) and then taking the large N 
limit. 
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Using the contour integral representation of Laguerre polynomials 

dw 



k T k / 2\ 

a L„,(—a 



aw k—l/-\ I —l\m+k 

- — re w [1 + aw ) 

2m 



(4.6) 



the summation over i in (4.5) can be easily performed. In this way, we arrive at the contour 
integral representation of the normal matrix integral (4.1) 






2Nw) 



N 



2Nw 



In this expression, it is easy to take the BMN limit 



J2 ^ 

N,J ^ oo g2 = —,\' = — fixed 
In this limit the integral (4.7) becomes a non-trivial function of g2 and A' 



(W^Oj) 



2 32^ f dw J JVV^U,,^,,,-!^ / 92/^ 

--e 8 



,^il^^i;^e^("+"'")|e^ 



1 



This integral is evaluated around the saddle point of the term proportional to J 



d„ 



{w + w ^) + logw 







w* 



'1 



1 1 



A' 



Then we find 



1 



{WaOj) = -Ae' 



where S and A are given by 



v^ 



A' 



A 



sinh 



1 + A'-i)-ief (1+^^)' 



1 \1 




f(l + ^ 


/I + A') 
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[i + VT+y 



(4.7) 



(4i 



(4.9) 



(4.10) 



(4.11) 



(4.12) 



Since (72 is the genus- counting parameter in the BMN limit, this expression gives an all-genus 
result of OPE between a single winding Wilson loop W^ and the 1/2 EPS BMN operator 
TtZ'^ . The expression of S agrees with the result of [37]. 

On the other hand, the coupling to the large winding loop W^"^^ {w ^ I) behaves very 
differently. As we saw above, the saddle point of z-integral is at z = z^ (2.15). Therefore, 
the OPE coefficient behaves as 



{W^^^Ok) ~ e 



2NKV^^TT+2{N+k) sinh^l k 



(4.13) 
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We can write down the OPE between the Wilson loop in the K^^ anti-symmetric repre- 
sentation and the chiral primary Oj in the form of the contour integral as before 

(Wa.Oj) = ™ ^ if[^e^^^ (z, + ^Vdet A (4.14) 

where K x K matrix A is given by 

( -y _l_ VX\N~J ( I V\\J _ yJ 

^jl — — 



1 Zi ~r 27Y ^1 

A, = lfi±^ivi3 i (,^1) (4.15) 

^« ^ 2N ] 

When K is of order 1 and A^ is sent to infinity, the dominant contribution comes from the 
diagonal elements of the matrix A 

K i^ _ / ./T\ ■^ X 



K rJidzi v2., f \^ 



= K{Waf'\WaOj) (4.16) 

Namely, the OPE coefficient factorizes when K is small. 

The more natural basis of chiral primary is given by the Schur polynomial of Z = 
diag(2;i,---,2;Ar) 

,. det(.^^^-^-) det(.^^^-^-) 

^"= det(.f-) = A(.) (^-^'^ 

Here hi > h2 ■ ■ ■ > h^ are the row length of the Young diagram R. Using the normal matrix 
model, one can easily show that the two-point function of Or is diagonal in R. Similarly, 
one can show that the OPE between W^ and the chiral primary Or is non-zero only if R is 
a single hook i.e., hi > 1, h2 = ■ ■ ■ = hk = I, hk+i = ■ ■ ■ = h^ = 







{WaOR)=0 if i?^y (4.18) 

Explicitly, the OPE is given by 



A^ / 2niz \ 2Nz^ 



{N + hi)'-^lMg^^{N + hi)) (4.19) 



N 

where k is the length of the first column and |i?| = /ii + /c — 1 is the number of boxes. It is 
interesting that A^ is effectively replaced hj N + hi. 
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For a more general Wilson loop Wji, the multi-hook chiral primary can couple to it. 
The computation becomes harder as the number of hooks increases. Therefore, we restrict 
ourselves to the case of single-hook chiral primary. The OPE between Wilson loop in the 
K^^ anti-symmetric representation and the single-hook chiral primary is easily obtained as 



dzi v\_,^ __-, / V\\ 



{Wa^O^) = / n ||e^" ^i' ^1 + ^ ^^^^ ^4-20) 



where K x K matrix B is given by 



Af-fc 



V^\ _ -N+k 



5., = lfi±^43 i (,^1) (4.21) 



2N/__'^J 
2N "'J 



Zi ~r oAr ^ 



Again, when K is of order 1 the dominant contribution comes from the diagonal part of B. 
In this case one finds 

{WA^Orrm) ~ K{Wa)''-'{WaOrrin) (4.22) 



As in the case of the Wilson loop correlator {Waj^), we expect that the OPE coefficient 
(Wai^Or) fails to factorize when K r^ N . We leave this analysis as an interesting future 
problem. 

Appendix C: Normal Matrix Representation of OPE between Wil- 
son Loop and Chiral Primary 

The OPE between Wilson loop and chiral primary is given by a complex matrix model 

/'[rf22]e-T'^(^^)^Tre5(^+^)C^(^) (4.23) 

Here, for notational simplicity, we have set 2^ = 1. This parameter can easily be restored 
by the rescaling of the matrix variable. Eq. (4.23) can also written as integrals over the real 
and imaginary parts of the complex matrix z, 

f[dXdY]e-^'^'''+^'^^Tie''OR{X - tY) (4.24) 

Then we can formally do the integral over the imaginary part Y to get 

|[rfXrfF]e-T'-(^'+^')i^Tre^e-^^'-^"^0^(X) 
= |[rfX]e-T'-^'lTre^e-^TK^)'o^(X) (4.25) 
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Note that we have taken into account the Wick contractions of powers of Y which had ap- 
peared in the chiral primary by operating with a derivative operator. Next we can diagonahze 
the Hermitian matrix X by conjugating it with a unitary matrix U 

X = f/diag(a;i,---,xjv)f/"^ (4.26) 

In the eigenvalue basis, the matrix Laplacian is given by 

where u; = U^^dU is the "angular" part. Since Oji{X) is independent of U, the Calogero-type 
term in the matrix Laplacian does not contribute. Therefore, the matrix integral reduces to 
the eigenvalue model 

l[dx,e-^'A{xf-Y.^^^irr^^'''^'-'''^i^)^Ri^) (4-28) 

for real eigenvalues Xi. Now, using the identity for an arbitrary function f{x) 

e~-^^"f{x) = e^^V (^^) e-^"' (4.29) 

the integral (4.28) becomes 

j\{dx^e-lT.,^'A(x)]^J:f On (i^) A (^^^ e'^E.^? (4.30) 

Now, we notice that the combination ^'^' = aj is a creation operator and e~2E»'^» is the 
wavefunction of the harmonic oscillator vacuum. The first factor and the last factor in(4.30) 
are both Slater determinants of the ground state wave-function of a set of A^ fermions moving 
in a harmonic oscillator potential, see eq. (2.19), up to a normalization factor 

A/'A(a:)e-^E."? = AT'A ( ^^ ] e~^E«"? = ^(xi, ■ ■ ■ , x^) (4.31) 

Using the bracket notation as in section 2.4, the integral is written as 

1^1 t 

(^lT7Ee^^""^""^^i?(«^)l^) (4-32) 

^^ j=i 

Finally, inserting the completeness relation for coherent states, (4.32) becomes a normal 

matrix model in the same way as we saw in section 2.4. 

This derivation of the normal matrix model shows that this relation is valid for any chiral 

primary Or. On the other hand, the Wilson loop part contains a normal ordering constant 

in general, hence the equivalence between complex matrix model and normal matrix model 

holds only for a certain representation where this normal ordering factor is common to all 

terms in the representation. For instance, the normal matrix model is exact for the OPE 

computation of the Wilson loop in the anti-symmetric representation. 
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5 Discussion 

In this paper, we have shown that the expectation values of circular Wilson loops and the 
OPE between Wilson loops and chiral primaries both can be written as normal matrix models. 
This representation as a normal matrix model is exact for some representations of Wilson 
loop and it is approximate for some other cases. 

This normal matrix model gives a unified description of 1/2-BPS chiral primary operators 
and the 1/2-BPS Wilson loops. In particular, there is a nice analogy between dual giant 
graviton and large winding Wilson loop. Namely, the dual giant graviton wrapped around 
S^ inside AdS^ is identified as a particle outside the droplet, whereas the large winding 
Wilson loop also corresponds to an isolated eigenvalue outside the droplet. Moreover it has 
been proposed that the object dual to a large winding loop is a D3-brane with world- volume 
AdS2 X S'^ and extended into AdS^. The direction that the dual object blows up into {i.e., 
AdS^ direction) is the same as the case of the giant graviton. 

In the case of giant graviton wrapped around an S^ inside S^, the droplet picture is a hole 
inside the droplet and the corresponding operator is a Schur polynomial of anti-symmetric 
representation. The corresponding Wilson loop is in an anti-symmetric representation of 
gauge group, and the dual object is a D5-brane with world-volume AdS2 x S"^. The S*^ 
part of the world-volume is inside S^, so the blowing-up direction {S^ in this case) is again 
common for both the chiral primary and the Wilson loop. However, the droplet picture of 
an anti-symmetric Wilson loop is a bit different from the giant graviton: it corresponds to 
a line dividing the droplet. Still there is a similarity between the droplet picture of giant 
graviton and anti-symmetric Wilson loop, namely they are both represented by objects inside 
the droplet. It would be nice to push this analogy further. 

We have noted that, when A is large, a symmetric representation Wilson loop is practically 
identical to a large winding number single trace Wilson loop. This explains why they can 
both be associated with the D3-brane. We also observe that, there is another phase which 
is stable for smaller A and a phase transition at A ~ 5.5. In that phase, the free energy 
resembles the D5-brane. It is interesting that this case has a D5-brane blowing up into the 
AdSs direction. 

We also considered the OPE between Wilson loop and chiral primary. At least technically, 
our normal matrix model gives an efficient way to compute this OPE. We found a peculiar 
selection rule for the OPE. For example, the Wilson loop in the fundamental representation 
couples only to the chiral primary with single hook representation. Similarly, the K^^ anti- 
symmetric loop couples to the chiral primary with at most K hooks. We should emphasize 
that this selection rule is exact. The meaning of selection rule in the bulk AdS side is not 
clear to us. It would be nice to understand this from the string theory side. 

Recently, the supergravity dual of the 1/2-BPS Wilson loop has been constructed [4, 10]. 
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The solution is characterized by a boundary value problem in two-dimensions, and the one- 
dimensional boundary line is divided into black and white regions. It is suggested [4] that 
these regions corresponds to the eigenvalue distribution of a hermitian matrix model. It 
would be nice to clarify this relation, if any. Our normal matrix model description seems 
to suggest the existence of M-theory version of the supergravity dual, where the boundary 
region is two-dimensional and the boundary data are droplets. 

We mainly considered the Wilson loop in the representation whose Young diagram has 
a number of boxes of order A^. We expect that when the number of boxes becomes of 
order A^^, the dual object can no longer be regarded as a probe of the fixed background 
^^5*5 X S^, rather the geometry is significantly deformed. In our droplet picture, the effect 
of such big Wilson loop will be to deform the shape of the circular droplet. For example, 
the anti-symmetric representation of the large winding loop (W^J) might correspond to 
two disconnected droplets when Kw ~ A^^. It would be interesting to find the eigenvalue 
distribution for this kind of large Wilson loop. 
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